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ABSTRACT 


We report the detection of a Cold Neptune rupianet = 21 ± 2 orbiting a 
0.38 Mq M dwarf lying 2.5-3.3 kpc toward the Galactic center as part of a cam¬ 
paign combining ground-based and Spitzer observations to measure the Galactic 
distribution of planets. This is the hrst time that the complex real-time protocols 
described bv lYee et al.l (j2015al) . which aim to maximize planet sensitivity while 
maintaining sample integrity, have been carried out in practice. Multiple survey 
and follow-up teams successfully combined their efforts within the framework of 
these protocols to detect this planet. This is the second planet in the Spitzer 
Galactic distribution sample. Both are in the near-to-mid disk and clearly not 
in the Galactic bulge. 


Subject headings: gravitational lensing: micro 


1. Introduction 

The 2015 Spitzer microlensing campaign is the hrst project with the specihc goal of 
characterizing the Galactic distribution of planets. Such characterization has three require¬ 
ments: 

1) A survey that is sensitive to planets in substantially different Galactic environments. 

2) Well-characterized selection. 

3) Ability to determine, at least statistically, the Galactic environment of each potential 
host, whether a planet is detected or not. 
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At present, microlensing is the only possible method to attack this problem becanse all 
other planet search techniques fail criterion (1). By contrast, microlensing is about equally 
sensitive to planet-hosting lenses at all distances along the line of sight from the Sun to the 
Galactic-bulge sources that are monitored for microlensing events. 


Initially microlensing planet searches were conducted in a fairly opportunistic way, dom¬ 
inated by targeted follow-up of “interesting” events. Because “interest ing” was not neces¬ 
saril y well-dehned , it wa s difficult to satisfy criterion (2). Nevertheless iGould et ahl (120lOl) 
and iGassan et al.l (120121 1 were able to construct subsamples with well-characterized selec¬ 
tion. In parallel, the RoboNet and MiNDSTEp t eams developed robotic a lgorithms flHorne 


Tsapras fc Snodgrass 20091: iDominik et al.ll2010l: iHnndertmark et al.ll2015l) to carry out the 


target selection in a repeatable and well-characterized manner. With the advent of second 
generation surv eys, particularly OGLE-IV and KMTNet, un iform selection is becoming rou¬ 


tine (see, e.g. ISuzuki et al.l 120151: IShvartzvald et ah 11201511 . In their high-cadence zones. 


these wide field surveys can obtain dense enough coverage that additional follow-up observa¬ 
tions are not necessary to detect and characterize planets. Hence, events can be monitored 
uniformly with a pre-defined observing sequence. 


The main problem has always been to determine the location of the lenses that are 
being probed for planets, including both those in which planets are detected and those for 
which they are not. For the great majority of microlensing events, the lens (host) star is not 
definitively identified and for many it is not detected at all. For these cases, the microlens 
parallax 


TTe = 


'^rel M 


9^ = kM TTreh 


AG 


K = 


~ 8 . 1 - 


mas 


( 1 ) 


0E /i’ """ c^AU Mq 

is the key to determining the lens mass and distance. Here, Tr^ei = AU(Zljj^ — and ^ are 
respectively the lens-source relative parallax and proper motio n, 9-p, is the angu l ar Ein stein 
radius, and M is the lens mass. As detailed by Figure 1 from iGould fc Hornel (120131) . the 
amplitude tte = vrrei/6*E is set by the fact that motion by the observer of 1 AU displaces the 
apparent angular separation of the lens and source by an angle vTi-ei, which in turn induces 
microlensing effects according to how large this displacement is compared to the Einstein 
radius. In principle, there are two ways to observe the displacement caused by parallax. 
First, the observer can wait to be moved by Earth’s orbital motion, creating a displacement 
relative to simple rectilinear motion between the source and lens. Second, two observers at 
substantially different locations can compare their observations. 


Some micro lens parallaxes have b een measured from the ground using the first effect. 


particularly for long timescale events (jPoindexter et al.l 12005) and including fortuitously a 


significant number of microlens planetary events (IGould et af 


20101) . However, the sample of 


events with parallaxes measured from the ground is extremely heavily biased toward nearby 
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lenses because they have projected Einstein radii te = AU/tte ~ 2-5 AU, roughly 3 times 
smaller than is typical of bulge lenses. Hence, this sample is almost unusable for measuring 
the Galactic distribution of planets 

The alternative is space-based parallaxes. At ~ 1 AU from Earth, Spitzer is ideally 
located to be such a “microlens parallax satellite”. It does, however, face a number of 
challenges. First, due to Sun-angle viewing restrictions, it can observe Galactic bulge targets 
(which are near the ecliptic) for only 38 continuous days out of the 8 months they are 
visible from Earth. Second, microlensing events must hrst be detected from the ground 
and identihed as reasonably planet sensitive before they are targeted by Spitzer, and thes e 
uploads occur 3-10 days before the observations begin (Figure 1 of lUdalski et ahl 120141) . 
Since microlensing events often evolve quite rapidly, these practical constraints intrinsically 
restrict the pool of targets. Finally, Spitzer observes at A = 3.6 /rm, roughly 4.5 times the 
wavelength of ground-based mi crolensing s earches . This creates technical problems, some of 
which are described below. See lYee et al.l (j2015al) for a more complete discussion. 


In 2014 the Director granted 100 hours for a pilot program to determine whether Spitzer 
could function e ffectively as a paral lax satellite. Special protocols were de veloped to rapidly 


upl oad targets Judalski et ah I2014J) . New techniques were developed b v lYee et ahl fl: 
and Calchi Novati et ah (2015) to break the famous “4-fold degeneracy” f Refsdai 19661: 


1994) 

that 

Goulc 

1995 


20141) . New techniques were developed b v lYee et ahl fl2015blj 


Gould 


l ad been believed to re quire extremely aggressive observational strategies (e.g 
Gandi fc Gould! 1 19971) . For the two events with measured (OGLE-2014-BLG - 


1050. IZhu et al.ll2015al) or strongly constrained (OGLE-2014-BLG -0124, Uda 
^E (using the standard technique for events with caustic features, IYoo et ah 


ski et al.ll2014l) 


20041) . the lens 


mass M and distance Dl were determined using inversions of Equation (fTj), 

0E „ „ AU 


M = 


TTrel — ^'eTTe; 


Dr = 


( 2 ) 


KTTe TTrel -|- 'nS 

Here is the source parallax which is almost always well understood because the source is 
visible. More critically, since the overwhelming majority of single-lens events do not show 
caustic features, a robust method wa s developed for estimating distances kinematically for 
events with measured tte but not ^e fiGalchi Novati et al.l 120151) . 


However, the 2014 Spitzer campaign was entirely focused on demonstrating the feasi¬ 
bility of the method, a nd no systematic ef fort was made to hnd planets. Nevertheless, one 
planet was discovered fjUdalski et al.ll2014l ). 


The successes of the 2014 campaign and the breakthroughs it precipitated led to the 
realization that satellite parallax measurements could be used to measure the Galactic distri¬ 
bution of planets. In 2015, 832 hours we re awarded for a program whose primary objective 
was to do just that fiGould et al.l 120141) . In fact, as specihcally argued in the proposal. 



































































several such annual campaigns will be required to acquire sufficient statistics to make this 
measurement. 


As discussed in detail bv lYee et ahl fj2015al) . the observational protocols required to 


1) maximize the sensitivity of the survey to planets, while 

2) maintaining well-characterized selection 


are both intricate and complex. We describe in some detail how these applied to the case of 
the planetary detection repor ted below . Howev er, those readers interested in a full under¬ 
standing must actually study lYee et ahl fj2015al) . 


Here, we report the discovery and characterization of the plane t OGLE-2015-BLG - 
0966Lb. This is the second planet (after OGLE-2014-BLG-0124Lb, Udalski et ah 2014 ) 
in the statistical sample of microlens-parallax planets that can be used to determine the 
Galactic distribution. The observations, and the observational protocols that guided them 
are described in Section |2l The ground-based and Spitzer lightcurves are combined to yield 
microlensing parameters for the event (including the microlens parallax tte) in Section |3l 
There are some challenges in the estimate of 6e relative to the usual case, which are discussed 
in Section 01 After resolving these, we present the system physical parameters in Section 01 
Finally, we discuss some implications of this discovery in Section [HI 


2. Observations 

2.1. OGLE Alert and Observations 


All 2015 Spitzer targets were chosen based on microlensing alerts issued by the Optical 
Gravitational Lens Experiment (OGLE) or the Microlensing Observations in Astrophysics 
(MOA) collaboration, with a substantial majority coming from OGLE. 


The OGLE alert for OGLE-2015-BLG-0966 was issued on 2015 May 11, well in advance 
of the Spitzer campaign, for which the hrst observations were on June 6. It lies at equatorial 
coordinates (17:55:01.02, —29:02:49.6), with Galactic coordinates (0.96, —1.82), placing it in 
OGLE held BLG505, which implies that it is observed at 20 minute c adence with OGLE ’s 
1.3m Warsaw Telescope at the Las Gampanas Observatory in Ghile (lUdalski et al.N2015l) . 
These dense observations permitted OGLE to alert the event, using its Early War ning System 
(EWS) real-time event detection software flUdalski et ahl Il994 lUdalskil 120031) . when the 
source had just entered the Einstein ring, and so was just 0.38 mag above its / = 19.62 
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baseline. OGLE continued to observe with this cadence throughout the event except for 
interruptions due to weather and monthly passage of the Moon through the bulge. 


2.2. Spitzer Observations 


A key element in meeting Criterion (2) for measuring the Galactic distribution of planets 
is that events must be selected for 5pzteer observations without allowing any knowle dge about 
the pre sence or absence of planets to influence that decision. As described in lYee et ah 
f)2015all . Spitzer targets can be chosen “objectively” or “subjectively”. If they meet certain 
specihed criteria as of 6 hours before target submission (Mondays at UT 15:00), then they 
are objective and must be chosen for observations with a certain specihed cadence. In this 
case, all the planets discovered (as well as all planet sensitivity, i.e., ability to detect planets) 
whether from before or after the Spitzer observations begin, can be incorporated into the 
analysis. In addition, if the event has not yet been selected objectively, the team may at any 
time still choose the event “subjectively” guided qualitatively by the goal of maximizing the 
sum (over all choices) of the products Yhi ^iPi where Si is the sensitivity of event i to planets 
and Pi is the probability that a microlens parallax will actually be measured. Only data 
taken (or rather, made public) after this selection date may be considered when calculating 
Si- The cadences of subjectively chosen events can also be specihed subjectively, but as a 
practical matter they are usually specihed to be the same as for objectively chosen events. 
Note that exclusion of Spitzer data from the analysis applies only to determining whether 
or not an event enters the sample. Once the event has a sufficiently well measured parallax 
to enter the sample, all Spitzer data can be applied to improve the precision of the parallax 
measurement or to search for planets. 


If there remains time for additional Spitzer observations after scheduling all targets ac¬ 
cording to their adopted cadences, this is applied to high-magnihcation events, with cadence 
ranked according to predicted magnihcation in the observation interval. 

A point of direct relevance to the present case is that if an event is initially chosen 
subjectively but later meets objective criteria, the objective status takes precedence. This 
means that all planets and planet-sensitivity from before the subjective alert can now be 
included in the analysis, but also that only observations following the objective alert can 
be used to determine if the parallax is measured well enough to enter the sample. If these 
objectively-based parallax data are not adequate for a measurement, then all the Spitzer 
data can be included, but then only planets and planet sensitivity from after the subjective 
alert can enter the analysis. 
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Figures 2 and 3 from lYee et al.l fj2015al) show flowcharts for this decision-making process. 


On Monday June 15, the Spitzer team chose OGLE-2014-BLG-0966 for “secret” ob¬ 
servations at 1-day cadence. The purpose of “secret” observations is to resolve the tension 
between two aspects of the observations. First, any event selected for Spitzer observations 
must be observed and continue to be observed at some predetermined cadence for a pre¬ 
determined amount of time. However, Spitzer targets can only be updated once per week. 
Hence, if the future course of an event is uncertain at the time the targets must be sent to 
Spitzer ior observations, it may be selected as “secret.” Until it is formally selected, none of 
the Spitzer observations may be included in the calculation of parallax but if the event turns 
out to have low planet sensitivity and/or poor chances of a viable parallax measurement, no 
commitment has been made to continue observing it in the future, which would be a waste 
of resources. Hence, “secret” observations permit the team to defer this choice until more 
information about the event is available, with a maximum “loss” of one week’s observations. 
In this case, however, the team gained sufficient confidence in the event just two days later 
to announce it publicly on 2015 June 15 UT 21:31. The cadence was specified as “objective”, 
meaning once per day, which could not then be altered. 


However, the next week, on Monday June 22, the event qualified for bonus observations 
because it was predicted to be high magnification (as seen from Earth) during that week of 
Spitzer ob s ervatio ns, June 24 UT 11:43 to July 1 UT 16:46. According to the prescriptions of 


Yee et ahl fj2015al) all observing time that is not allocated for regular observations should be 


applied to high-magnification events, rank ordered by the 1 a lower limits of their predicted 
peak magnification during the observing interval. This prediction is for magnification at 
Earth if that is all the information available (as was the case for this event in that week). 
OGLE-2015-BLG-0966 was estimated as A^ax > H and so was assigned 8-per-day cadence. 
It received a total of 58 observations during that week. In fact, this week did prove to contain 
the event’s A^ax ~ 100 peak as seen from the ground. The magnification of the event during 
the next week was less than the cutoff for distributing extra Spitzer observations, so it was 
observed at the previously determined cadence of 1 /d ay. This week it a lso met the objective 


criteria for selection of rising events (criteria “B” in lYee et al.ll2015al) . This meant that if 


the event parallax could be adequately measured using Spitzer data from after July 1 UT 
16:46, then all planets and planet sensitivity could be included. Otherwise only planets 
and sensitivity from after 2015 June 15 UT 21:31 could be included (assuming all Spitzer 
data were enough to measure a microlens parallax). For the final “week” (actually about 
10 days), the event again met the criteria for additional high-magnification (as seen from 
Earth) observations and was slated for a cadence of 4-per-day. It may seem strange that 
it would fail these criteria on the fifth week but pass them on the sixth week, since the 
event was falling during this time. However, the criteria for receiving extra observations 
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due to high-magnification are (necessarily) based on predicted magnihcation, not actual 
magnihcation, and these predictions improve with time. Second, the competition from other 
events varies each week, so there is strict correspondence between (predicted) magnihcation 
and observations only within a given week, not between weeks. 

Altogether, Spitzer observed this event a total of 129 times, each with 6 dithered 30 s 
exposures. 


2.3. MOA data 

MOA independently identihed this event on 16 June and monitored it as MOA-2015- 
BLG-281 using its 1.8m telescope with 2.2 deg^ held at Mt. John New Zealand. In contrast 
to most other observatories, which observe in I band, MOA observes in a broad R-I band 
pass. The MOA cadence for this held is 15 minutes. Thus, during the long June-July nights, 
OGLE and MOA together nominally cover this held for about 18 hours per night. However, 
the weather in New Zealand is far worse than in Chile, so the fraction of nights that actually 
have this near-continuous coverage is under 50%. 


2.4. Ground-Based Follow-up Data 

2.4-1- Follow-up Strategy 


As discussed in lYee et al.l (l2015al) , ground-based follow-up strategy is intimately con¬ 
nected with event selection, both objective and subjective. Approximately 150 (nominally) 
point-lens events were selected for S'pifzer observations during the 2015 campaign, which far 
exceeds the resources of all follow-up groups combined to densely monitor events to search 
for planets. This tension has two interrelated implications. First, follow-up groups are ex¬ 
plicitly encouraged not to monitor events that are heavily monitored by surveys, with survey 
coverage being described in some detail at the time of announcement of each event. Second, 
there is a strong bias for the Spitzer team to select events that are heavily monitored by 
surveys, exactly because coverage is not dependent on limited follow-up resources. Indeed, 
for events (such as OGLE-2015-BLG-0966) that have 20-minute OGLE cadence, there is an 
extremely strong bias because additional follow-up from the same time zone would be redun¬ 
dant, and hence the follow-up observing resources are best applied to other events without 
such intense survey coverage. Similarly for events in the 16 deg^ core KMTNet helds, which 
have roughly 15-minute cadence from three observatories, these considerations apply even 
more strongly. OGLE-2015-BLG-0966 lies approximately in the middle of one of the four 
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KMTNet prime fields. However, the Spitzer team’s map of these helds was precise enough 
to recognize that the event lay in a gap between chips, so that it was not actually covered 
by KMTNet. (KMTNet did not at this time publicly list the events that it was monitoring). 


As OGLE-2015-BLG-0966 approached its peak (HJD' 7205.2, July 1.7), the very high ca¬ 
dence and high quality of OGLE data permitted an accurate estimate of the peak magnihca- 
tion Ajnax ~ 90, which would make the event highly sensitive to planets (IGriest fc Sahzadeh 


19981 ). Because OGLE would normally densely cover the event during the long Ghile night. 


while MOA would cover the even longer New Zealand night, the need for follow-up would 
have appeared minimal. This is particularly true for the fairly large number of follow-up 
telescopes in Ghile, whose observations would be completely redundant with OGLE. 


However, the peak of this event happened to occur when the Moon was passing through 
the bulge, during which time OGLE does not observe. Hence, several follow-up groups 
concentrated their efforts on this event. Nevertheless, the event did not gain the undivided 
attention of follow-up groups due to several competing events. Most notable was OGLE- 
2015-BLG-0961, which peaked less than 12 hours earlier and at even higher magnihcation 
^max = 200. However, this event was covered by KMTNet and so assumed much lower (but 
still not zero) priority from follow-up groups. 

In particular, we note that the planet was discovered only because follow-up groups 
recognized that there would be no survey coverage of this event over Ghile. 


2 . 4 . 2 . LCOGT 


Las Gumbres Observatory Global Telescope Network (LGOGT) provided ground-based 
observations primarily from its southern ring of 8 1.0m telescopes s ited at GTIO/Ghile, 
SAAO/South Africa and Siding Spring/Australia flBrown et al.N2013[) . Two telescopes in 
Ghile are equipped with the new generation of Sinistro imagers that incorporate 4kx4k 
Fairchild GGD-486 B1 GGDs and offer a held of view of 27' x 27'. All other 1.0m telescopes 
support SBIG STX-16803 cameras with Kodak KAF-16803 front illuminated 4096x4096pix 
GGDs, used in bin 2x2 mode with a held of view of 15.8' x 15.8'. Observations were also 
made from the 2.0m Faulkes North Telescope in Haleakala, Hawaii, using its lO'x 10' Spectral 
camera (Fairchild GGD-486 B1 GGD). All telescopes in the network are equipped with a 
consistent set of hlters. SDSS-f' was used for the large majority of these observations, with 
a small number made using the Bessell-K hlter. Due to the constraints described abo ve, 
LGOGT employed its TArget Prioritization (TAP) algorithm flHnndertmark et al.ll2015[) to 
select a sub-set of events from the Spitzer target list based on their predicted sensitivity to 
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planets, which were drawn from Spitzer targets that fell in regions of lower snrvey observing 
cadence. Since OGLE-2015-BLG-0966 falls in a region of high snrvey cadence, it was not 
selected for observation nntil the Moon’s passage throngh the Bnlge interrnpted snrvey 
observations. At that point the event was flagged for high density observations on both 
1.0m and 2.0m networks. These observations were condncted as gronps of 2-10 exposnres 
repeated at intervals of a few minntes to honrs, with less dense observations being taken 
as the event retnrned to baseline. Hence, the data typically have dense packets of coverage 
followed by short gaps. 


2.4-3. Danish Telescope 


The Danish 1.54m telescope is one of the national telescopes hosted by ESO at La Silla 
in Ghile. After a snccessfnl refnrbishment in 2012 by the Gzech company Projectsoft, it was 
eqnipped wit h the first rontinely operated mnlti-color instrnment providing Lncky Imaging 
photometry (jSkottfelt et al.l 1201511 . The instrnment itself consists of two Andor iXon+ 897 


EMGGDs and two dichroic mirrors splitting the signal into red and visnal bandpasses. For 
the 2015 microlensing campaign the camera was oper ated at 10 Hz, and ln cky exposnres 
were calibrated and tip-tilt corrected as described by iHarpscie et ahl (12012h. Photom etry 


was obtained from the collapsed images nsing the DanDIA pipeline flBramich II2008I) and 
based on rontines of the RoboNet rednction pipeline. A modified version of ARTEMiS 


flDominik et al.l 120081) was deployed to coordinate the observation of microlensing targets, 
and in the case of OGLE-2015-BLG-0966 resnlted in confirmation of the anomaly from a 
second site. 


2 . 4 . 4 . CTIO-SMARTS 

For the dnration of the Spitzer campaign, the Microlensing Follow Up Network (/iFUN) 
donbled its normal allocation to 6 honrs per night on the dnal optical/IR ANDIGAM camera 
on the 1.3m SMARTS telescope at GTIO. This facility was tasked with several objectives that 
were not always completely compatible. These inclnded regnlar monitoring of Spitzer targets 
in low-cadence OGLE fields in order to predict their fntnre behavior, sparse monitoring of 
all targets in order to measnre their if-band sonrce flnxes nsing ANDIGAM’s IR channel, 
and dense monitoring of events that were at fairly high or very high magnihcation in order to 
detect planets. Hence, similar to the LGOGT telescopes, GTIO-SMARTS observed in blocks 
with short gaps dnring which other goals (inclnding other science and shared resonrces) were 
pnrsned. 
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2.5. Other Bands 

OGLE, LCOGT, and GTIO-SMARTS all obtained G-band observations in order to 
measnre the V — I sonrce color in order to characterize the sonrce. As described above, 
GTIO-SMARTS automatically obtained if-band observations through the IR channel simul¬ 
taneously with both the R-band and J-band optical observations. These also are used for 
source characterization and are not included in the analysis. 


2.6. Data Reduction 


All ground-based data were reduced using stan dard algorithms. All data entering the 


maiii analysis used variants of image s ubtraction flAlarc 
2013 ). GTIO-SMARTS used DoPhot f Schechter et al 
characterization analysis. 


fc Luptonl Il998l: iBramich et al 


19931) reductions for the source- 


As explained in lYee et al.l (l2015al) . no previously existing Spitzer reduction software 
was suitable for variable stars in crowded helds. Therefore customized software had to be 
developed. This will be described in a forthcoming work (Galchi Novati et al. 2015, in prep). 


3. Lightcurve Analysis 


Real-time lightcurve analysis, both manual and automated flBozza et al.ll2010[l was con¬ 
ducted during the event, and immediately alerted observers to the presence of the anomaly. 
Extensive offline analysis was subsequently carried out once data gathering on the event was 
complete. 


The analysis reported below is based on a simultaneous £t to ground-based and Spitzer 
data. However, the key event characteristics can be most easily understood by considering 
the two lightcurves separately, and indeed many characteristics can be inferred by visual 
inspection. 


Ignoring for the moment the post-peak perturbation centered at fpert = 7205.75, the 
ground based lightcurve shown in Figure [T] peaks at to = 7205.19 almost 4.85 mag above 
baseline, indicating a high-magnification event A^ax ^ 85 (i.e., more if the baseline source is 
blended). In the neighborhood of the peak, such events are characterized by flux evolution 
F{t) = Fpeak(l + (^ — where G ff = uotK, tp is the Einstein timescale and uq 

is the impact parameter normalized to 6 -^ flGouldl 119961) . Inspection of the peak region 
shows feff = 0.68 days. (For example, GTIO data at ±1.50 days from peak are 0.96 mag 
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below peak. Hence, tes = 1.5 day(10°-®*°'®® — 1)“^/^ = 0.68 day.) The fact that there is 
no pronounced dip just prior to the abrupt rise due to a caustic crossing at tec = 7205.64 
shows that the companion/host mass ratio is very small, i.e., this is a planetary rather 
than binary system. Hence, the center of mass of the system is very close to the “center of 
magnihcati on”. It is the sou rce motion relative to the “center of magnification” that produces 


the overall iPaczvhskil (1l986l) curve over the unperturbed portions of the light curve. Hence, 
the caustic that induces the perturbation, and which lies along the planet-host axis, points 
toward this center of magnification. This implies that the source is moving at an angle 
a = tan“^(tefr/(tpert — io)) = 0.88 radians (50.5°) relative to the planet-host axis. There 
is no dip between the caustic entrance and exit, which implies that the caustic edges are 
separated by signihcantly less than the source radius p (normalized to ^e). Hence, the 
source radius crossing time, = pts, can be estimated ~ (tec — ^pert) *sin(Q;) = 0.07 days. 
While it is not obvious by inspection, a simple point-lens £t to the unperturbed lightcurve 
shows that the source is essentially unblended, so that in fact Hmax = 85, so uq = 0.012, 
and hence = tefr/wo = 57 days, and thus p = 1.2 x 10“^. That is, (to, mo, ^e, o, p) = 
(7205.19, 0.012, 57day, 51°, 1.2 x 10“^). These by-eye estimates agree reasonably well with 
the parameter values derived from the best-fitting models (described below), presented in 
Tables [U and El The remaining two parameters that can be extracted from the ground-based 
lightcurve, i.e., the planet-star mass ratio q and the planet-star separation s (normalized to 
^e), cannot be estimated by eye and must be determined from the fit. 

The effect of adding the Spitzer data is to determine the microlens parallax tte, which 
is incorporated directly into the fit in equatorial (north, east) coordinates. However, as 
illustrated in Figure 1 of iGouldl fjl994l ). the impact of these measurements is better visualized 
in the frame defined by the projected Spitzer-Earth axis D^, in which 

fo,© to,sat _ 


ATI 

tte = — (At, A/1); 


At = 


in 


A/1 = ±Mo,® - ±Mo, 


sat ? 


( 3 ) 


and where the subscripts indicate parameters as measured from Earth and Spitzer. 


For the intuitive analysis of the Spitzer lightcurve (Fig. [T]) we begin with the external 
information that the timescale Ie — 57 days is known from the ground-based lightcurve. In 
fact, since Earth-Sptteer motion (projected on the sky) is about 20kms“^ at the time that 
we will evaluate Equation (jS]), the timescales will not be exactly the same. However, this 
is expected to be a relatively small effect. Hence, we ignore it here with the proviso that 
we will later check for self-consistency. The Spitzer lightcurve shows a factor ~ 4 rise over 
the 30 days of observations ending at 7222.14, where it is clearly turning over, implying a 
peak fo.sat ~ 7225. Thus, the interval between peaks, normalized to the Einstein timescale 
is At = (to,® — to,sat)/^E ~ 0.35. Since 30 days before the peak, Wsat ~ 0.5 {A ~ 2), the 
magniheation at peak is at least Amax,sat ^ 8 (more if the Spitzer source flux is blended), so 
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I A/?| = 1^0,0 — Mo,sat| — Mo,sat ^ 0.12. SpitzcT was 1.39 AU from Earth midway between the 
peaks, with = 1.27 AU after taking account of projection. Since it is approximately due 
West of Earth, the two coordinates in Equation ([3]) basically correspond to East and North, 
respectively. Hence, we derive tte.e ~ —0.28, and |vrE,Ar| < 0.1 As a hnal check, we note that 
these values imply a projected velocity in the Earth frame of h = AU/(7rEtE) ~ 100kms“^, 
meaning that Earth’s motion is not completely negligible. In principle, we could make 
recursive corrections, but the main point of this exercise is to show that most of the results 
of the full light curve analysis can be basically derived by simple inspection of the Spitzer 
and ground-based lightcurves. 


We systematically model the combined Earth-based and Spitzer lightcurves, notwith¬ 
standing the above analysis showing that hve of the seven Earth-based parameters are well- 
determined by visual inspection and simple analysis. We conduct a systematic search of 
parameter space using two different techniques, one based on an (s, q, a) grid and the other 
based on lightcurve morphologies. Models are evaluated based on the reduced of fh® fif 
per dataset {Xred)- 'f'o account for variations in the estimation of photome tric errors in dif¬ 
ferent datasets, we adopt the common practise (e.g. 
the errors according to: 


Bachelet et ah (1201211 ) of re-normalize 


^new 


ao 



(4) 


Coefficients qq and ai are set such that the Xred of each dataset relative to the model 
equal unity. 


We found that the only viable solutions do in fact have the hve parameters as approx¬ 
imately specihed above. Moreover, s and q are well localized, except that (as is often the 
case) there is a degenerate solution s flGriest fc Sahzadehlll998l) . 


Figures [2] and [3] illustrate the geometries of all eight solutions, and Table [T] gives the 
parameters for these solutions. The notation (-(-, —) is used to indicate solutions where lens- 
source relative trajectory approaches the caustic from the positive or negative 6 direction 
in the lens-plane geometry. Two indices are given for each solution: one for Earth-based 
observations and one for Spitzer. According to Table [H the best wide and close solutions 
differ by < 1. Hence, this very common degeneracy is completely unbroken in the 
present case. However, since the microlensing parameters of these two sets of solutions are 
very similar, their physical implications are basically the same. 


Within each set of solutions (wide or close) the physical implications are even closer, to 
the point of being nearly identical. Therefore, we present in Table [2] a simplihed summary 
of each set of solutions, with the parameter values being averages over the four solutions 
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and the “errors” taking acconnt of both the £t errors at different minima and the differences 
between minima. 


4. Einstein Radius Estimate 


We estimate 6 *e nsing the standard approach flYoo et ahl |2004[) . i.e., estimating the 
sonrce snrface brightn ess from its deredden ed color {V — I) q and the calibrated color/surface- 
brightness relation of iKervella et ahl (120041 1. and then comparing this to dereddened magni- 
tnde Isfi- 

We begin (as is usnal) by assnming that the sonrce is behind the same dust column as 
the bulge clump stars. The field shows substantial differential reddening, so we restrict to 
a 60" radius about the source, wherein the clump appears compact. Using OGLE data, we 
End that the source lies AJ = 3.15 mag below the clump, whereas using CTIO data we find 
AJ = 3.01. This difference is typical of the difhculty in centroiding the clump in the vertical 
(magnitude) direction. We adopt AJ = 3 .08 T 0.07 and so Igfl = /ci,o + = 17.54 where 

/ci,o = 14.41 is adopted from Nataf et ah (2013). 


Determining {V — I)s (and so (U —/)s,o) presents greater challenges given that the source 
is quite faint at baseline V ~ 21.9 and that the full moon was passing through the bulge 
when the source was most highly magnified. OGLE systematically acquires U-band data 
for all its fields, but of course avoids observing near the Moon. Hence, the OGLE (V — J)o 
color estimate has relatively large errors: A(U — /)* = —0.29 ± 0.06. GTIO deliberately 
targeted the event near peak (and so Moon passage) in order to obtain high signal-to-noise 
measurements. Some of these observations were corrupted by scattered moonlight, but most 
appear fine. These data yield A{V — I)s = —0.35±0.03. We combine these OGLE and GTIO 
measurements by standard error weighting and obt ain (V — I),^ a = (V — I)ci,o + A(U — /) = 


0.72 where we have adopted (V — I)ci,o = 1-06 from lBensbv et ahl (1201311 . If we assume that 
the source is exactly at the distance to the clump and adopt M/= —0.12, then this implies 
[M/, (V — /)o]s = (2.96, 0.72). This is a plausible pair of values for a star just entering the 
sub-giant branch from the turnoff. That is, the assumption that the source lies behind all the 
dust is self-consistent, since it implies that the source inhabits a reasonably well-populated 
part of the color-magnitude diagram. 


However, another interpretation is that the source lies on the upper main sequence (since 
its color is similar to the Sun). Then, Mj ~ 4.15, so that the star lies 1.2 mag in front of the 
Bulge in distance modulus, or at Ds ~ 4.5 kpc from us. In this case, it would still lie about 
125 pc below the Galactic plane and so behind most but possibly not all the dust toward 
the Galactic bulge. 


We first estimate the source radius under the assumption that it is a bluish bulge 
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subgiant and then address how this estimate may be a ffected if the source is in fact a disk 
main sequence star. After using iBessell &: BrettI fjl988l) to convert (V — I) to (/ — K), the 
above procedure leads to estimates, 


Q 

6 ^ = 1.07 ± 0.10 /ras; = — = 0.76 mas; 

P 


Pge 


Q 

= — = 4.8masyr“b (5) 


If the source were a disk main sequence star, but nevertheless lay behind all the dust, 
then these estimates would not be affected at all. Only the dereddened color and magnitude 
enter the calculation, and these would not change. If the source does he in front of some of the 
dust, then it is both intrinsically redder in [V — I) and fainter in I than the above calculations 
would imply. Of course fainter stars have smaller radii at hxed surface brightness, but redder 
stars have lower surface brightness (so larger radii) at hxed magnitude. Hence, these two 
effects tend to cancel. Since the total amount of dust behind the source is unlikely to be 
large and the two effects tend to cancel, and since at this point we cannot reliably estimate 
how much dust does lie behind the source, we will proceed on the assumption that the source 
is behind all of the dust. However, in Sectional we discuss how this issue could be partly or 
fully resolved in the future. 

5. Physical Parameters 


The extraction of physical parameters is made simpler by the “resolution” (or rather 
irrelevance) of the ambiguity of the amplitude of A(3 (Section |3]), but more complicated by 
the ambiguity in the source distance (Section 0]). We begin therefore with the lens mass 
estimate, which does not depend on the source distance 

Ov 

M = -^ = 0.38 ± 0.04 Mq] rupianet = qM = 21 ± 2M® (6) 

KTTe 

Similarly, for the projected velocity in the geocentric and heliocentric frames 

Vgeo(A^, E) = ^ ^ = (0, -124)kms"^ Wei = Vgeo + v®,± = (0, -95)kms“^ (7) 

where v®^j_(A^, E) = (—0.8, 28.8)kms“^ is Earth’s velocity projected on the plane of the sky 
at the peak of the event and where we have ignored the slight differences among the four 
wide solutions. 

Regardless of the source location, the relative parallax is well-determined 

TTrei = = 0.19 mas, (8) 

but depending on the distance to the source, this leads to two different distance estimates 
for the lens Dl = 3.3kpc (bulge source) or = 2.5kpc (disk source). However, within 






19 


the context of our program of determining the Galactic distribution of planets, these two 
outcomes are basically similar: foreground disk lens. In Section [6l we discuss how this 
ambiguity may be resolved with future observations. 

Finally, the ambiguity 'm. Dl gives rise to an ambiguity of the same size in the star- 
planet projected separation = sO-^Dl = 2.7AU (bulge source) or 2.1 AU (disk source). 
Adopting a snowline r^naw = 2.7 AU {M/Mq), these separations are at r_i_/rsnow = 2.2 and 
1.7 snow-line distances, respectively. Hence, in either case, this planet is a “cold Neptune”. 
Note that in making these estimates, we have adopted s = 1, i.e., midway between (and 10% 
different from) the close and wide solutions. 

6. Discussion 


This is the second microlensing planet (after OGLE-2014-BLG-0124Lb), whose mass 
and distance have been characterized with the aid of parallax observations using Spitzer. 
However, it is the hrst to be so characterized as part of a program specihcally designed 
to measure the G alactic distribution of planets by means of well-defined selection criteria 
flYee et al.ll2015al) . While naively such well-defined criteria might seem to lead to clear cut 
observing decisions, the example of OGLE-2015-BLG-0966 demonstrates that they in fact 
create a set of complex interlocking constraints, both hard and soft, on decision processes 
of multiple semi-autonomous decision makers working toward a common goal. Specifically, 
there was the struggle to determine whether surveys would cover this event, including both 
static but non-obvious (KMT) information and dynamic OGLE information, impact on bal¬ 
ancing with coverage of other events (also impacted by survey coverage), and wild oscillations 
in Spitzer cadence due to objective criteria. See Section [21 


In the present case, these follow-up observations that were subject to these interlocking 
constraints led to detection and characterization of a planet while preserving the integrity of 
the selection process. In subsequent papers we will show that similar sensitivity to planets 
(while maintaining sample integrity) was achieved under similar conditions for several other 
high-magnification events. 


The character of the planet, continues to confirm that “cool Neptunes are common” 
fiGould et al.ll2006l) . even though that original claim was based on just two detections. How¬ 
ever, the central focus of the present effort is not the frequency of planets as a function of 
mass, but of Galactocentric radius. Nothing can be said so far about this in part because 
there are so far only two planets in the sample, but mainly because the sensitivity of the 
surveys to planets as a func tion of Galactocentric radius has not yet been determined (al¬ 
though see lZhu et al.ll2015bl for initial work in that direction). Nevertheless, it is intriguing 
to note that both detected planets are in the near-to-mid disk. 
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While the lens is clearly in the disk, it is nnknown at this point whether the sonrce is in 
the bnlge or the disk. While the sonrce distance is of overall secondary interest, it does affect 
the lens distance and planet-star separation at the 25% level and therefore shonld 
be resolved if possible. We note that the heliocentric projected velocity Wei does tend to 
favor a disk sonrce jnst becanse the solntions all have the heliocentric projected velocities 
pointing basically ont of the Galactic plane, and indeed somewh at retrograde. By contrast, 
as expected naively (and confirmed bv ICalchi Novati et al.ll20151) events with bnlge sonrces 
and disk lenses tend to have Vhei aligned with Galactic rotation. This is becanse the Snn 
and the lens both basically partake of this rotation while bnlge sonrce proper motions tend 
to be both isotropically distribnted and of low amplitnde. 

The situation w onld be greatly clarihed by measuring the proper motion of the source, 
as Yee et ah (2015b) did for the Spitzer event OGLE-2014-BLG-0939. In contrast to that 
case, however, in which the source was bright and easily distinguished from neighbors, which 
enabled precise astrometric measurements from over a decade of OGLE observations, the 
OGLE-2015-BLG-0966 source star is both too faint and too close (0.7") to a neighbor for re¬ 
liable astrometry. However, using difference image astrometry applied to high-magnification 
images, the current position of the lens is measured with a precision of about 20 mas. Hence, 
a single epoch of high-resolution imaging 10 years from now should be able to reliably de¬ 
termine whether the source proper motion is consistent with disk motion (typically about 
6 mas yr“^). 
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Table 1. All solutions. 


Parameters 

(+, +),wide 

(+, —),wide 

(—, +),wide 

(—, —),wide 

(+, +),close 

(+, —),close 

(—, +),close 

(—, —),close 

X^/dof 

15136/15278 

15138/15278 

15137/15278 

15135/15278 

15136/15278 

15151/15278 

15154/15278 

15135/15278 

to 

7205.1979 

7205.1953 

7205.1955 

7205.1937 

7205.1978 

7205.1945 

7205.1936 

7205.1938 


0.0012 

0.0013 

0.0013 

0.0012 

0.0012 

0.0012 

0.0012 

0.0012 

UO 

0.01136 

0.01136 

-0.01157 

-0.01152 

0.01136 

0.01136 

-0.01162 

-0.01149 


0.00010 

0.00010 

0.00010 

0.00010 

0.00010 

0.00010 

0.00010 

0.00009 

tfi (days) 

57.7 

57.8 

57.7 

57.7 

57.7 

57.7 

57.5 

57.8 


0.4 

0.4 

0.4 

0.4 

0.4 

0.5 

0.4 

0.4 

p (10-b 

14.08 

14.00 

14.06 

14.06 

14.07 

14.02 

14.12 

14.04 


0.14 

0.15 

0.14 

0.13 

0.15 

0.15 

0.15 

0.14 

'TE,N 

0.0234 

-0.0561 

0.0397 

-0.0412 

0.0234 

-0.0538 

0.0317 

-0.0412 


0.0012 

0.0015 

0.0012 

0.0013 

0.0012 

0.0059 

0.0094 

0.0012 

'TE,E 

-0.238 

-0.237 

-0.232 

-0.237 

-0.239 

-0.264 

-0.272 

-0.237 


0.006 

0.012 

0.003 

0.006 

0.006 

0.007 

0.007 

0.006 

a (deg) 

50.7 

50.6 

-50.4 

-50.4 

50.6 

50.6 

-50.3 

-50.4 


0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

s 

1.1141 

1.1146 

1.1148 

1.1147 

0.9093 

0.9091 

0.9088 

0.9091 


0.0033 

0.0032 

0.0033 

0.0033 

0.0026 

0.0025 

0.0027 

0.0026 

q (10-4) 

1.67 

1.67 

1.68 

1.68 

1.70 

1.70 

1.71 

1.70 


0.05 

0.05 

0.05 

0.05 

0.05 

0.05 

0.06 

0.05 

Blend 

-0.027 

-0.024 

-0.027 

-0.026 

-0.027 

-0.025 

-0.031 

-0.024 


0.008 

0.009 

0.008 

0.008 

0.008 

0.008 

0.008 

0.008 


Note. 


HJD-2450000. 
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Fig. 1.— OGLE-2015-BLG-0966 lightcurve. Gombined data from 10 telescopes (color coded) 
trace the ground-based light curve nearly continuously. With the except ion of the ~_6hr 
post peak bump, it is well characterized by a high-magnihcation point-lens Paczyhsl^ (1986) 
curve. As described in the text, hve of the seven parameters (to, Mq, p) needed to 
describe this curve can be read off lightcurve or extracted with very simple analysis. The 
remaining two (s, q) (planet-star separation and mass ratio) require more detailed modeling. 
The Spitzer lightcurve is aligned to the OGLE scale so that eqnal “magnitndes” represent 
equal magnihcations. The microlens parallax tte can be well-estimated simply by comparing 
the Spitzer and gronnd-based light curves. See Section [3l (OGLE and MOA data are binned 
for display in the Figure, but not in the ht. Ground-based data points with nncertainties 
> 0.2 mag are snppressed in the hgure to avoid clutter, but are included in the ht.) 
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Fig. 2.— Thumb-nail sketches of all eight caustic geometries for planetary microlensing 
event OGLE-2015-BLG-0966. The black and red lines represent the trajectories as seen 
from Earth and Spitzer, respectively. The right ends of the red lines represent the end of 
Spzfeer observations at HJD 7222.14. All eight solutions have essentially the same physical 
implications. First, while the overall caustic structures (closed black curves) differ between 
close (left) and wide (right) binary solutions, the central caustics, which are the only part 
probed by the Earth and Fpifoer observations, are virtually identical. These are related by the 
s y-)- degeneracy, and since s ~ 1.1 is close to unity, they correspond to physically similar 
systems. The microlens parallax tte is essentially determined by the offset between Spitzer 
and Earth trajectories at the same epoch (magenta dashed line segments). See Equation ([3]). 
The amplitude of tte does differ slightly between the (±, ±) and (±, =f) solutions because 
the source is on the same side of the lens as seen from Earth and Spitzer for the former 
and on the opposite side for the latter, leading to a different distance perpendicular to the 
direction of motion. However, because the event is high magnihcation as seen from Earth 
(■*^ 0,0 ^ 1); this difference is itself small, and because the separation along the direction 
of motion is much larger than the separation perpendicular, this has almost no effect on 
the magnitude of the parallax, tte, which is what enters the main physical parameters. See 
Table [TJ The only substantial difference among these eight solutions is that the motion is 
somewhat north of west for the (-I-, ±) solutions and somewhat south of west for the (—, ±) 
solutions. However, this small difference has no impact on the main physical parameters, 
which do not depend on this direction. 
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ejB^ (East) 


Fig. 3.— Detail of caustic structure for one of the eight degenerate solutions shown in 
Figure [2] (close The colored circles represent the data points on the light curve 

shown in Figured] The black hlled circle on the trajectory represents the size of the source. 
As described in Section |3l the angle between source trajectory and the binary axis (same 
as the caustic axis) is determined directly by the time offset of the anomalous peak from 
the overall peak, relative to the curvature of the peak, without reference to any model. 
However, the derivation of (s, g), which determines the overall caustic structure in this plot, 
does require detailed modeling. See Section [3l 



Table 2: Abbreviated table of parameter values averaged over the four solutions. 


Parameters 

Wide solution 

Close solution 

Best 

15135 

15135 

y.a 

7205.196 ±0.002 

7205.195 ±0.002 

\uo\ 

0.0115 ±0.0001 

0.0115 ±0.0001 

tE (days) 

57.7 ±0.4 

57.7 ±0.4 

P (10-^) 

14.04 ±0.15 

14.07 ±0.15 

kE,N 

0.04 ±0.02 

0.03 ±0.02 

TJ'E.E 

-0.24 ±0.01 

-0.25 ±0.02 

|a| (deg) 

50.5 ±0.2 

50.5 ±0.2 

s 

1.115 ±0.004 

0.909 ±0.003 

o 

1 

1.68 ±0.05 

1.70 ±0.06 
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